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ABSTRACT 

We present a statistical study of the non-thermal X-ray emission of 27 young rotation powered pulsars 
(RPPs) and 24 pulsar wind nebulae (PWNe) by using the Chandra and the XMM-Newton observations, 
which with the high spatial resolutions enable us to spatially resolve pulsars from their surrounding 
PWNe. We obtain the X-ray luminosities and spectra separately for RPPs and PWNe, and then in- 
vestigate their distribution and relation to each other as well as the relation with the pulsar rotational 
parameters. In the pair-correlation analysis we find that: (1) the X-ray (2-10 keV) luminosities of both 
pulsar and PWN (Lx,psr and £x.pwn) display a strong correlation with pulsar spin down power E and 
characteristic age r, and the scalings resulting from a simple linear fit to the data are -Lx,psr oc ijo.92±o.o4 



and Lx.pwn oc E 



1.45±0.08 



(68% confidence level), respectively, however, both the fits are not statistically 



acceptable; (2) £x.psr also shows a possible weak correlation with pulsar period P and period derivative 
P, whereas Lx.pwn manifests a similar weak correlation with P only; (3) The PWN photon index Fp^n 
is positively correlated with Lx,pwn and Lx,pwn/E. We also found that the PWN X-ray luminosity is 
typically 1 to 10 times larger than that from the underlying pulsar, and the PWN photon indices span 
a range of 1.5 < Fpwn ^ 2. The statistic study of PWN spectral properties supports the particle wind 
model in which the X-ray emitting electrons are accelerated by the termination shock of the wind. 

Subject headings: radiation mechanisms: non-thermal - star: neutron - stars: pulsar: general - X-rays: 
general 



1. INTRODUCTION 

The rotation-powered pulsars (RPPs) are known as 
rapidly spinning and strongly magnetized neutron stars 
that are radiating at the expense of their rotational en- 
ergy. The X-ray emission of RPPs may contain both ther- 
mal and non-thermal components. The thermal emission 
might be further divided into non-pulsed and pulsed com- 
ponents. The non-pulsed component, originates from the 
cooling of the neutron star, is from the whole pulsar sur- 
face with a characteristic temperature of about 0.1 keV, 
while the pulsed component comes from hot spots 1.0 
keV) on the pulsar surface, which are heated by the bom- 
bardment of relativistic particles streaming back to the 
surface from the pulsar magnetosphere (Kundt & Schaaf 
1993, ZavHn et al. 1995, Gil & Krawczyk 1996). The 
non-thermal pulsar emission is from the pulsar magneto- 
sphere, and it might also contain pulsed (e.g., Cheng & 
Zhang 1999; Zhang & Harding 2000) and non-pulsed com- 
ponents (e.g., Tennant ct al. 2001; Becker et al. 2004). 
In some cases, a pulsar wind nebula (PWN) is found to 
surround a RPP. The X-ray emission of the PWN is non- 
pulsed and often dominates the non-thermal emission of 
the system. 

A lot of efforts have been devoted to the statistical 
studies of pulsar X-ray emission properties, with partic- 
ular emphasis on the efficiency of conversion of the pul- 
sar spin down power E into X-ray luminosity. By us- 
ing the Einstein data, Seward & Wang (1988) found that 
Lx oc E^'^^, where Lx is the 0.2-4.0 keV X-ray luminosity 
of the pulsar (plus PWN). Becker & Triimper (1997) ob- 



tained Lx — 10 ^E using a sample of 27 pulsars observed 
by ROSAT, where Lx is the total X-ray (0.1-2.4 keV) lu- 
minosity of the pulsar plus PWN. However, in these two 
works, the thermal emission may contribute significantly 
to the total X-ray luminosity given the adopted soft X- 
ray band. Saito (1998) analyzed 16 RPPs observed by 
ASCA (2-10 keV) and found Lx cx E^/^. Possenti et al. 
(2002) reported Lx oc E^-^'^ using 39 pulsars observed by 
ASCA, RXTE, BeppoSAX, Chandra and XMM-Newton. 
The X-ray luminosities in Saito (1998) and Possenti et al. 
(2002) also include the total emission due to the pulsars 
plus PWNe, given the limited spatial resolutions of ASCA, 
RXTE, and BeppoSAX. Cheng et al. (2004) divided the 
total X-ray emission into pulsed and non-pulsed compo- 
nents, and found that the X-ray luminosity of the pulsed 
one follows Lx.pui oc E^-'^, which agrees with the model 
prediction Lx oc E^-'^^ by Cheng & Zhang (1999). For 
non-pulsed emission, they got Lx,npui oc L'^ '*, where they 
supposed that the non-pulsed component comes mainly 
from PWNe and the contribution of non-pulsed compo- 
nent from pulsars is negligible. It is worth noting that the 
scatter in the relation is large, as pointed out by Possenti 
et al. (2002), who performed study including the estimates 
of the observational errors and showing that the linear fit 
is statistically unacceptable. 

All these previous works suffer from the low spatial 
resolution of the detectors, making it difficult to resolve 
the emission of the pulsars from that of their surround- 
ing PWNe. It is the purpose of our current work to re- 
solve and to analyze the pulsar and the PWN emission 
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separately. Thanks to the high spatial resolution obser- 
vations performed with Chandra and XMM-Newton, we 
have been able to satisfactorily investigate 27 pulsars and 
24 PWNe, for which we have determined the non-thermal 
X-ray fluxes and spectra in the 2-10 keV band. Then we 
have carried out separated statistic studies of RPPs and 
PWNe, and tested the consistence of their emission prop- 
erties with current models. The organization of this paper 
is as following: the sample and the data processing are 
presented in section 2; the statistical analyses of the X-ray 
spectral properties of RPPs and PWNe are given in sec- 
tion 3; we discuss the physical implications of our results 
in section 4 and summarize our work in section 5. 

2. SAMPLE AND DATA PROCESSING 

We collect pulsar and PWN samples from the obser- 
vations by Chandra and XMM-Newton, which both have 
high spatial resolutions, i.e., ^ 1" and ~ 6", respectively. 
We take the Chandra data directly from the literatures, 
and if there are no published results, we analyzed the data 
in this paper. The XMM-Newton data are adapted only 
if there are no relevant data from Chandra for the same 
source. All the XMM-Newton results are taken from litera- 
tures. Totally we obtain the X-ray spectra of 27 RPPs and 
24 PWNe. In our samples, millisecond pulsars (MSPs) are 
not included. It is generally believed that MSPs have ever 
undergone an accretion-driven spin-up phase and they are 
usually old and regarded as a significantly different class. 
Similar study on the MSPs is also limited by the rare data 
available. Therefore we do not analyze MSPs here, al- 
though we discuss them when compare our analysis with 
the previous work including MSPs. 

In our samples, there are 15, out of 27, spectra of pul- 
sars obtained from the archived Chandra data. We select 
only the pulsars detected by the Advanced CCD Imag- 
ing Spectrometer (ACIS) in the Timed Exposure (TE) 
Mode, in which a pulsar is able to be resolved spatially 
from its surrounding PWN. We calibrate the data with 
CIAO (ver 3.4) and CALDB (ver 3.3.0). We first repro- 
cess the Level 1 data for the correction of the charge trans- 
fer inefficiency (CTI) effects, then clean the background 
and remove the afterglow. Time intervals with anomalous 
background rates associated with particle flare events are 
further rejected in the Level 2 data. And then the pulsar 
positions are obtained by the celldetect tool in CIAO. Fi- 
nally, the spectra are extracted from the Level 2 data and 
then fit with XSPEC. We use both the power-law (PL) 
and the power-law-|-blackbody (PL-I-BB) models to fit the 
pulsar spectra. If the resulted spectral indices are consis- 
tent within errors in both models, then the results from 
the PL model are used, otherwise those from the PL+BB 
model are used. In our spectral analysis, we show errors 
at the 90% confidence level. 

Pileup occurs when more than one photon are collected 
in one pixel within a CCD readout frame, since those pho- 
tons can only be recorded as a single photon event whose 
energy is the sum of the collected photons. Therefore 
pileup may affect the results of spectral analysis. Accord- 
ing to section 6.14.2 in the Chandra Proposers' Observa- 
tory Guide v.lQ, the effect of pileup can be omitted if the 

* http:/ /cxc. harvard.edu/proposer/POG/ 
^ http:/ /cxc. harvard.edu/toolkit/pimms.jsp 



pileup fraction is <10%. However, pileup does affect the 
spectral analysis even if the pileup fraction is <10%. For 
example, since the pileup fraction of PSR J1930-I-1852 is 
estimated to be only 6%, its spectral index is reported to 
be 1.09^0 09 without pileup correction (Lu et al. 2002), 
whileas the spectral index is 1.35tQ;5o after pileup correc- 
tion (Camilo et al. 2002). In our spectral analysis, we 
first estimate the pileup fraction using PIMM^, and then 
add a pileup model in the spectral fitting of the pulsars 
if the pileup fraction is higher than 3%. Totally, there 
are 8 pulsars in which the pileup model is included in 
the spectral fitting, i.e., PSRs J0205-f6449, J0537-6910, 
B0540-69 (and its PWN), B0833-45 (Vela), J1747-2958, 
J1846-0258, J1930-hl852 and B1951-H32. 

The absorption column density (Nn) is obtained in 
several ways. (1) For 6 pulsars (PSRs J0205+6449, 
J0537-6910, J1747-2958, J1846-0258, J1930-I-1852 and 
B1951-I-32) with bright PWNe, Nn are obtained from 
the spectral fitting of their PWNe, and then fixed when 
fitting the spectra of the pulsars. (2) PSRs B0540-69 
and Jl 124-5916 are embedded in SNRs 0540-69.3 and 
G292.0-f 1.8, respectively, and their PWN spectra below 
2.5 keV are strongly affected by the SNRs. At the same 
time, constraining TVh with emission above 2.5 keV is dif- 
ficult because of the small absorption in this high energy 
range. Therefore their iVn are obtained by fitting the 
pulsar spectra that the contamination of the SNR emis- 
sion is negligible, and then are fixed in the spectral 
fitting of their PWNe. (3) The PWNe associated with 
PSRs B0355+54, J1617-5055, B1823-13, B1929-I-10 and 
J2229-I-6114 are not bright, and iVn is determined by 
jointly fitting the spectra of both the pulsar and its PWN. 
(4) PSRs J0633-t-1746 and B0833-45 have been studied 
extensively, and their iVjj values used in our spectral fit- 
ting are taken from Caraveo et al. (2004) and Pavlov et 
al. (2001). 

The luminosity uncertainty is crucial in our analysis of 
correlations, and should be considered carefully. Since the 
X-ray luminosity is given by Lx — Aird^ fx, where d is 
the pulsar distance and /x is the 2-10 keV X-ray flux, the 
Lx uncertainty should be derived from the uncertainties of 
both fx and d. The uncertainty of fx is derived from those 
of the normalization and the photon index in the spectral 
fitting. For the fluxes taken from literatures, their uncer- 
tainties are extrapolated from the published ones by the 
ratios of the fluxes in 2-10 keV to those in the correspond- 
ing published energy ranges. 

The distances are usually not well constrained, thus the 
distance uncertainty may dominate the luminosity uncer- 
tainty. There are several cases in our samples: (1) the 
distances of 7 pulsars are derived from the radio disper- 
sion measures, and their errors are conservatively taken 
to be 40%, as estimated by Cordes & Lazio (2001); (2) 
the distances of 14 pulsars are obtained via their associ- 
ated SNRs, and some of them are shown with published 
distance errors in literatures, while for the others without 
published errors a conservative error of 50% is taken; (3) 
PSRs J0537-6910 and B0540-69 are both located in the 
Large Magellanic Cloud (LMC), whose distance is taken 
as 50 kpc, and an error of 10 kpc for a conservative esti- 
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mation is adapted (Bradley 2007). 

We summarize the properties of all the 27 RPPs and 24 
PWNe in Tables 1 and 2, respectively. In Table 2 only 22 
PWNe are associated with the RPPs listed in Table 1. The 
other two pulsars are excluded from Table 1: (1) Camilo 
et al. (2004) suggested that PSR J1016-5857 is too faint 
to be resolved from its background PWN in the Chan- 
dra Observation; (2) Hesscls et al. (2004) found that the 
spectrum of PSR J2021+3651 can be fit by a BB model 
and is thus dominated by thermal components. On the 
other hand, 5 out of the 27 RPP samples in Table 1 are 
not listed in Table 2 for PWNe, because the following 5 
pulsars have no PWN reported: B0628-28, B0656+14, 
B0823+26, B0950+08 and B1055-52 (Ogelman & Te- 
pedelenlioglu 2004, Becker et al. 2004, De Luca et al. 
2005). 

3. ANALYSES AND RESULTS 

The pulsar photon indices are distributed in a range of 
1 ^ Tpsi. < 3 (as shown in Fig[T]). It should be noted 
that a significant fraction, about ~ 15% (4 out of 27), of 
the sources have soft spectra of Fpsr > 2, which may raise 
problems for current models as discussed later in §4. The 
photon indices of the PWNe span a narrower range (Fig 
[TJ. As discussed below (§4), this is consistent with the 
pulsar wind model. 

We investigate below the correlations between the X-ray 
emission properties of RPPs and PWNe, and between the 
emission properties and the pulsar rotational parameters. 
The rotation parameters include the period P, the period 
derivatives P, and some derived parameters, e.g., the mag- 
netic field B = 3.3 X 10^^{PP)^/'^G, the characteristic age 
T = P/2P, and the spin down power E = iir^IP/P'^, 
where a typical moment of / = lO'^^g cni^ is assumed. We 
have taken the values of P and P from the pulsar catalog 
by Manchester et al. (20050. 

In order to evaluate the significance level of the corre- 
lations of the two parameters concerned, we calculate the 
widely used Pearson correlation coefficient (r) , the Spear- 
man rank correlation coefficient (rg), and the two-sided 
significance level (ps) of the Spearman rank test. The re- 
sults are listed in Table 3. 

In addition to the correlation tests, we also perform a 
linear fit using the least square method (LSM) to the rel- 
evant relations of the parameter pairs. Since the fitting 
results are usually dominated by a few data points with 
much smaller observational errors than the others, we also 
perform a linear fit without the observational errors for 
comparison. The fitting results are all listed in Table 3, 
and shown in the relevant figures. In the following we 
present the results in details. 

3.1. Correlations between the RPP emission properties 
and the pulsar rotational parameters 

We study the RPP emission first. Strong correlations 
appear between the X-ray luminosities of pulsars (ix,psr) 
and the pulsar rotational parameters (see Table 3 and Figs 
[5] and [3]). First, Lx,psr is negatively correlated with r and 
positively correlated with E, which are supported by the 
Spearman tests: rg = —0.81 and Ps < 0.0001 between 



Lx.psr and r; and r^ = 0.82 and Ps < 0.0001 between 
-^x,psr and E. We also note that there are some hints 
of the correlation hold for Lx,psr vs P and P separately, 
with the relevant Spearman rank correlation coefficients 
of —0.66 and 0.69 respectively and both significance levels 
< 0.001. The correlations between Lx,psr vs P and P will 
disappear when the sample includes both the MSPs and 
the normal RPPs, just as shown by Possenti et al. (2002). 

Despite the Pearson and Spearman correlation coeffi- 
cient support the existence of a correlation between Lx.psr 
and E (or r), a simple linear fit to the logarithm of the 
data points with the observational errors included does not 
produce a statistically acceptable model. In fact, it results 
(here Lx,psr and E are in units of erg and r in years; 
see also Figs 2 and 3). 

Lx,psr = 10-0-8±1-3£;0-92±0.04(^2 ^ 3.6) 
Lx,psr = 10381±O V-119±" 05(;^2 ^ 4 9) 

Here and elsewhere in this paper, the uncertainties on the 
linear fits are reported at 68% confidence level. Previous 
authors (notably Possenti et al. 2002) also noticed that 
a large scatter in the plot prevents to obtain an accept- 
able fit of the data with a simple power law dependence 
of Lx,psr on E. Hence this relation must only be seen as 
an empirical average trend and not suitable for predicting 
the luminosity of any specific source. 

We have also explored a linear fit which does not account 
for the uncertainties on the values of Lx,psi- It turns out 

T* 1 n-4.2±3.7 T^l.OiO.l 

T* _ 1 n38.9±0.9_-1.4±0.2 
-'^X.psr — T 

(see Figs [5] and [3]) . A comparison between the current 
ix,psr versus E relation and the previous studies is shown 
in Fig [21 It can be seen that the relation we obtain above 
is close to the one between the pulsed X-ray emission and 
the spin down power in Cheng et al. (2004), which indi- 
cates that most of the non-thermal X-ray emission from a 
pulsar is pulsed. 

As already done by Possenti et al. (2002) using a sam- 
ple including also the MSPs (but not disentangling PWN 
from RPP emission), we also try to fit Lx with aP^P'^. 
This gives the relation 

ix,psr = (40± l)p-3.4±0.3p0.77±0.07(^2 ^ 3.5). 

The nominal result of this (still statistically unacceptable) 
fit would suggest a preferred dependence of I/x,psr on E/P: 
however we note that, accounting for the uncertainties on 
the parameters, the simpler dependence on E (recovered 
in the work of Possenti et al. 2002) is also viable. 

We also study the pulsar spectral properties, and check 
if there is any correlation between the pulsar spectral in- 
dex Fpsr and the pulsar rotational parameters. Inspection 
of Fig. [31 may indicate the occurrence of a positive corre- 
lation of Fpsr with P and r and of a negative correlation 
of Fpsr with P and E. However, a numerical test indicates 
that all these correlations are too weak (the Spearman co- 
efficients |rs| are all < 0.60, see Table 3) for drawing any 
firm conclusion with the available data. 
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3.2. Correlations between the PWN emission properties 
and the pulsar rotational parameters 

We study here the correlations between the PWN X-ray 
himinosity Lx.pwn and the pulsar rotational parameters. 
As shown in Table 3 and Fig. O a weak positive correla- 
tion between Lx.pwn and P has been detected. Table 3 
and Figs [2] and [5] also show that Lx.pwn is strongly corre- 
lated with E and r. The linear fits with the observational 
errors result 

Lx^p^„ = 10-19.6±3.0^1.45±0.08(^2 ^ 3.7) 

Again, ix.pwn and E are in units of erg s~^. Like for 
the emission from RPPs, the adopted linear model in the 
logarithm of the data does not provide a statistically ac- 
ceptable fit. 

Trying a linear fit which does not account for the uncer- 
tainties on the values of ix.pwn, we obtain 

T* _ in-14.9±6.0 pl.3±0.2 

-'^X.pwn — -L'J ^ 

T* — in40.5±l.l^-1.7±0.3 

''^X,pwn — -LU ~ 

We note that the slope of the relation Lx.pwn °^ ^^'^ 
somewhat different from that of the pulsar, p^^. cx E. 
The same holds true comparing Lx.pwn oc E^'^^ with 
^x.psr oc E'^-^^. It is worth noting that, as seen in Fig. 
m the scatterings in the relation of Lx.psr versus E and 
that of the PWNe are comparable and both are large. 

As seen in Fig [5] and Table 3, there is no evidence for 
strong correlations between Fpwn and the pulsar rotational 
parameters. The Spearman rank test also supports this 
eye-ball study. However in Figs [1] and [7| we see obvious 
positive correlations between the photon index Fpwn and 
the X-ray luminosity Lx,pwn or the X-ray conversion efh- 
ciency Lx,pwn/E. We will discuss the physical implication 
in §4. 

3.3. Correlations between non-thermal emission 
properties of RPPs and PWNe 

For those samples with both the RPP and the PWN 
non-thermal X-ray emission measured, we test the corre- 
lations between them. As shown in Fig[8l a strong corre- 
lation between the X-ray luminosity of RPPs and that of 
PWNe appears in our samples, while no correlation shown 
between their photon indices. 

The correlation test between two luminosities gives rg — 
0.91 and Ps < 10"'', and trying a linear fit to the relation 
leads to Lx,pwn = lO^^'^'^^-^ix.p"'^- Given the strong 
positive correlations of Lx.psr and ijx,pwn versus E sepa- 
rately, a strong correlation between the two luminosities 
might be naturally expected. However, the slope of the 
relation between the two luminosities is somewhat differ- 
ent from that expected from the previous two relations, 
ix.psr E and ix.pwn oc E^-^. This may be because 
that the samples used are different. For example, those 
data points with E < lO'^'^ergs s"'^ are not included for 
the relation of Lx,pwn versus E since no obvious PWNe 
were detected, and they seem to result in a smaller slope 
for Lx,psi versus E relation (Fig [J). However, it should 



be noted that both Lx,psr and Lx,pwn are actually modu- 
lated by the source distance. The above correlation might 
be due to the effect of distance modulation. 

In our samples, the X-ray luminosity ratio between 
PWNe and RPPs, as shown by /x,pwn//x,psr in Fig. ^ 
varies in the range of 0.1 — 30, about 2 orders of magni- 
tude, and PWNe are generally brighter than their related 
RPPs, typically ix,pwn/-Z^x,psr 1 — 10. Fig. 9 also tells 
us that a more energetic pulsar does not tend to transfer a 
bigger fraction of E into PWN emission than to the pulsar 
emission, and vice versa. 

Gotthelf (2003) reported a linear relation between Fpsr 
and Fpwn for the Crab- like pulsars, but we find no correla- 
tion between them in our samples (Fig. [5]). We also check 
the relation using the samples of Gotthelf (2003), and a 
similar relation appears. These might suggest that the lin- 
ear relation exists only in the Crab-like pulsars, which are 
very young. 

4. DISCUSSIONS 
4.1. Non-thermal X-ray luminosities of RPPs and PWNe 

The detected X-ray emission from the RPPs and their 
PWNe is powered by the pulsar rotation energy. In our 
sample, the conversion efficiency of E to the 2-10 keV X- 
ray emission varies in the range of (Lx.psr + Lx.pwn)/E 
10~^ — lO""', with the mean at ~ 10"'^, so usually only a 
small fraction of the spin-down power goes into the non- 
thermal X-ray emission. The non-thermal X-ray emission 
in a source is usually dominated by the PWN rather than 
the pulsar, and on average, (ix,pwn/-^x,psr) ~ 10. This 
implies that the relations of the pulsar X-ray luminosity 
and the spin-down power obtained in the previous works 
using the low spatial resolution observations at > 2keV 
might be dominated by the PWN emission. 

In this work we can separately analyze the luminosity 
and the spin-down power relations for the pulsars and their 
PWNe, thanks to the high spatial resolution of the obser- 
vations. A strong positive correlation between Lx.psr and 
E is obtained in this paper, similar to the previous results, 
as shown in Fig [2l We compare our results with the pre- 
vious work in the following. However, one should keep in 
mind that there are obvious differences in the analysis pro- 
cesses, i.e., we can separate the RPP and the PWN emis- 
sion and do not include the MSPs in the sample, while the 
previous work did not separate the RPPs and the PWNe 
and included MSPs in the analysis. 

The relations for the X-ray luminosity of the RPPs in 
2-10 keV band are Lx.psr oc ii;0-92±0.04 (uncertainties on 
ix.psr included) and -Z^x psr °^ £;i o±o.i ^^q^ accounting 
for the uncertainties on ix.psr), respectively. They are 
roughly in agreement with the scaling found by Becker & 
Triimper (1997), who used the X-ray luminosity in the 
0.1-2.4 keV band. The Lx — E relations obtained by 
Saito(1998) and Possenti et al. (2002) appear steeper than 
both our derived relations for RPPs, but they are roughly 
consistent with both Lx,pwn oc iji -isio.os ^^^^^ L'x p^^ oc 
j^i.3±o.2 (Fig [2]). This may suggest that their relations 
could also be influenced by the PWN emission due to the 
lower spatial resolution. The relations Lx,pui oc £;i-2±o.08 
and ix,npui oc i;i-4±o.i^ obtained by Cheng et al. (2004) 
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with ASCA data, are similar to our results. 

We have shown that the weak negative correlation of 
^x.psr versus P and the weak positive correlation of Lx,psr 
versus P lead to a strong positive correlation between 
-^x.psr and E in our MSP-excluding RPP samples. On 
the other hand in the previous work including the MSPs 
and the normal RPPs, the correlations between the X- 
ray luminosity (might include the PWN emission) and P 
or P disappear whereas a trend of Lx versus E is still 
there (e.g., Possenti et al 2002). Moreover, the MSP sam- 
ples alone also obey a similar correlation (Possenti et al. 
2002). All these factors together strongly suggest that the 
X-ray luminosities of RPPs, including the MSPs, are only 
dependent on their spin-down powers. 

Although there is a strong correlation between Lx.psr 
and E, the scattering in this relation is large and the lin- 
ear fit with the observational errors included usually gives 
a statistically unacceptable result, as suggested by Pos- 
senti et al. (2002). Lx,psr at given E may spread over 2-4 
orders of magnitude, as seen in FiglH The uncertainty in 
the distance determination and the momenta of inertia are 
not expected to lead to such a large span, so other intrinsic 
factors may work, e.g., the viewing angle effect, etc. The 
scattering in the ix.pwn — E relation is comparably large, 
which is somewhat strange, since the PWN emission is less 
influenced by the viewing angles. 

4.2. Non-thermal X-ray spectra of the RPPs 

There are mainly two scenarios to produce the non- 
thermal X-rays in the magnetospheres of pulsars. The 
outer gap scenario (e.g., Cheng et al. 1998; Wang et 
al. 1998; Cheng & Zhang 1999) produces a downward 
synchrotron-curvature cascade, where the secondary elec- 
trons/positions produce X-rays by synchrotron emission. 
Another scenario is the polar gap scenario, e.g., Zhang 
& Harding (2000) proposed the "full polar cap cascade", 
where the non-thermal X-rays are produced by resonant 
inverse Compton (IC) scattering off the thermal X-ray 
photons. In both scenarios the Lx.psr oc E relation is gen- 
erally predicted, although the X-ray spectra are not easy 
to understand. 

Wc note that a significant fraction of pulsars with very 
soft spectral indices, Tpsr ~ 2 — 3, may pose questions 
on current models. In a cascade, the monoenergetic pri- 
mary electrons emit monoenergetic curvature photons, 
which subsequently turn into still monoenergetic pairs 
in a soft photon bath. The fast energy loss of the sec- 
ondary pairs in the magnetic field produces synchrotron 
emission, which have an photon spectrum with a power 
law index Fi = 1.5. If the cascade continues the pho- 
tons produce next-generation pairs and then synchrotron 
photons with index r2 = 1 -I- ri/2 = 1.75; furthermore. 
Pa = 1 -I- r2/2 = 1.875... So the indices will never be big- 
ger than 2, in contrast with the soft spectra. Actually this 
discussion could also work if IC rather than synchrotron 
emission is involved since the index of synchrotron and IC 
emission is equal for the same energy distribution of pairs. 
As for the polar gap scenarios, the synchrotron emission 
at X-rays is weak because the secondary pairs with small 
pitch angles produce synchrotron emission well above the 
cyclotron frequency in the strong pulsar magnetic field. 



typically ~ 100 keV. Zhang & Harding (2000) proposed 
that the X-ray emission is dominated by the low energy 
tail in the resonant IC emission. However this tail may 
be hard with index P < 2, as shown in some Monte Carlo 
simulations (e.g.. Fang & Zhang 2006), although the cases 
might be more complicated when more factors such as the 
viewing angle are taken into account. 

Wang & Zhao (2004) reported the possible negative cor- 
relations between Fpsr and f2 and between Fpsr and Ci 
where ^ is the generation order parameter characterizing a 
pulsar under the scheme of cascade processes (Zhao et al. 
1989; Lu et al. 1994; Wei et al. 1997). A similar negative 
correlation between Fpsr and ( in anomalous X-ray pulsars 
and softer gamma-ray repeaters had also been reported 
(Marsden & White 2001; Lu et al. 2003), suggesting that 
a common mechanism may operate in both normal and 
anomalous pulsars. These observational results seem in 
contrast with the predicted positive correlation between 
Fpsr and C by Lu et al. (1994). Here, we also check the 
relation between Fpsr and the generation order parameter 

= 1 + (0.6 - (ll/14)logP + (2/7)logA5)/l-3 (Eq. 6 of 
Wang & Zhao 2004) and list the correlation results in Ta- 
ble 3. It turns out that although there may be some hints 
of such a negative correlation in our sample, the correla- 
tion tests do not support it strongly, Ts = —0.5. So the 
current data are not good enough to test the theoretical 
predictions of Lu et al. (1994). 

4.3. X-ray spectra of PWNe 

In the standard Kennel & Coroniti (1984ab) model for 
the Crab nebula, the young Crab pulsar loses its rotational 
energy predominantly in the form of a highly relativis- 
tic particle wind, which encounters with the surrounding 
medium in a termination shock and become visible by syn- 
chrotron emission downstream from the shock. In this con- 
text the energy in the relativistic wind is transferred into 
post shock magnetic field and accelerated particles with 
energy distribution of Ne{Ee) oc E~p. Chevaher (2000) 
discussed the PWN spectra with emphasis on the cooling 
of the X-ray emitting electrons, which leads to a steeper 
index p-\-l for high energy and fast-cooling electrons, and 
hence a spectral transition of synchrotron photons from 
(p -\- l)/2 in the slow-cooling regime to {p 2)/2 in the 
fast-cooling regime at break frequency (vc)- 

The data fit with a single power law model to the indeed 
broken power law would result in a spectral index always 
in the range of (p + l)/2 < Fpwn < (p + 2)/2. An electron 
index value p w 2.2 is generally obtained in theoretical 
works on particle acceleration in relativistic coUisionless 
shocks, by both numerical calculations (e.g., Achterberg et 
al. 2001) and analytic analysis (e.g., Keshet & Waxman 
2005), and also inferred from observation in other kinds 
of astrophysical relativistic shocks, e.g., GRB afterglows 
(e.g., Freedman & Waxman 2001). Our results show an 
narrow index range of 1.5 < Fpwn ^2.1 unless one source 
with somewhat higher value ~ 2.5, suggesting an electron 
index of p ~ 2.2. This consistence with particle wind mod- 
els gives a strong support to the Fermi-shock acceleration 
in PWNe. 

We show that the PWN spectral parameters are not 
strongly correlated with the pulsar rotational parameters 
(Fig[S]). Gotthelf (2003) reported the correlation between 
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Fpsr and Fpwn for nine Crab- like pulsars. Our studies show 
that such a correlation is probably not a common prop- 
erty for all RPPs. Therefore, the electron spectrum and 
its evolution in a PWN are not determined by the central 
pulsar, consistent with wind models where the emission 
comes from a relativistic shock between wind and envi- 
ronment interaction. 

The relation of Fpwn with PWN luminosity Lx.pwn and 
the conversion efficiency Lx,pwn/E (Figs [1] and [7] and Ta- 
ble 3) could be understood qualitatively in the framework 
of pulsar wind models taking into account the electron 
cooling effect on spectral profile (e.g.. Chevalier 2000). 
If pulsar loses most of its rotation energy through par- 
ticle winds, then higher E corresponds to stronger cooling 
and hence lower spectral break j/c, which further means a 
larger index Fp^n in a fixed observational energy range. 
In the meantime, a higher E corresponds to a larger 
^x.pwn, no matter i/c is below or above the observational 
range, and corresponds to constant X-ray conversion effi- 
ciency for fast cooling regime {vc below observed range) 
or larger Lx.pwn/E for slow cooling regime (i^c above ob- 
served range). Therefore we have softer PWN spectra 
(larger Fp^n) for more luminous PWNe (larger Lx,pwn) 
and higher energy conversion efficiency {Lx,pwn/ E). This 
consistence supports the wind-shock model for PWNe. In 
this context, the transition of Fp„,i from high to low values 
in Fig. [T]suggests that the spectral break Vc locates at 2-10 
keV for Lx,pwn ~ lO^'^ergss"-^. This may give constraint 
to wind model parameters. 

5. CONCLUSIONS 

In this work, using the available samples of 27 RPPs 
and 24 PWNe observed by Chandra and XMM-Newton, 
we obtain the non-thermal X-ray spectral properties, i.e., 
luminosities and spectral indices, of RPPs and PWNe sep- 
arately. We then analyze their distribution and correlation 
with each other and with pulsar rotational parameters. 

• As to the correlations we find: (1) Lx^psr and 
^x.pwn display a strong correlation with both E 
and t; (2) Lx.psr also shows a possible weaker cor- 
relation with P and P, whereas ix.pwn manifests 
a similar weak correlation with P only; (3) Fpwn is 
positively correlated with ix.pwn and the efficiency 
of conversion of rotational energy loss in X-ray lu- 
minosity Lx,pwn/E. 



• Trying to fit the logarithm of the data with a sim- 
ple linear fit, we find: Lx psr = 10-0-8±i-3£;0.92±o.04 

and ix.pwn = 10-19-6±30i;1.45±0.08^ jjow- 

ever, both the fits are statistically unacceptable. 
Not accounting for the uncertainties on the ob- 
served luminosity, the aforementioned relations be- 
come L^^p,, = I0-4.2±3.7^i.o±o.i L^^^^^ = 
2o-i4.9±6.0£;i.3±o.2^ respectively. Since the scatter 
in the relation for PWN (whose emission should be 
less affected by viewing angle) is comparably larger 
than that for RPPs, the scatter in the relation is 
more probably intrinsic to the sources. 

• The PWN X-ray luminosity is typically 1 to 10 
times larger than that from the underlying pulsar. 

• The pulsar photon index spans a range of 1 < 
Tpsr ^3. A significant fraction of RPPs with low 
E show soft spectra of Fpsr > 2, which seems not 
consistent with the current models and urges for 
further investigation of the non-thermal X-ray emis- 
sion mechanisms of pulsars. 

• The PWN spectral properties are consistent with 
the particle wind model: the photon index range 
1-5 ^ Fpwn ^ 2 is consistent with that expected 
from the shock- accelerated electrons of index p ~ 2; 
the correlations of Fpwn with Lx,pwn and the con- 
version efficiency Lx,pwn/E are consistent with the 
wind model; no correlation between Fpwn and the 
pulsar rotational parameters also implies that the 
cooling process is not related to the center pulsars 
but to the interaction of the pulsar wind with its 
environment. 
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Table 1 

Rotational parameters and emission properties of 27 pulsars 



PSR Name P P TVh d Tp^, /x (2-lOkeV) L^^psr Det. Ref. 

s ss^i lO^^cm^^ kpc ergss~icni~2 ^-i 



J0205+6449 


0. 


,066 


1, 


,939 


X 


10- 


-L3 





B0355+54 


0. 


,156 


4, 


,397 


X 


10- 


-15 





B0531+21 


0. 


,033 


4. 


.228 


X 


10" 


13 





J0537-6910 


0. 


,016 


5, 


,178 


X 


10- 


■14 





B0540-69 


0. 


,050 


4, 


,791 


X 


10- 


13 






B0628~28 


1, 


.244 


7, 


,123 


X 


10- 


■15 


0, 


1 q+0.03 

■-■^"^-0.02 


1, 


Al\ 


J0633+1746 


0, 


.237 


1, 


.097 


X 


10- 


■14 


0, 


,001 


0, 


,25; 


B0656+14 


0, 


.385 


5, 


.500 


X 


10- 


■14 


0. 


Q40+0.002 
■"^'-'-0.002 


0, 


,29" 


B0823+26 


0, 


.531 


1, 


,709 


X 


10- 


■15 


0, 


r, + 0.088 

■"-0.0 


0, 


3+' 


B0833-45 


0, 


.089 


1, 


.250 


X 


10- 


-13 


0, 


,017 


0, 


,29" 



0l42 3:2^^^^ TF^^ 7.4+^-^ X 10-13 9.1+''^Y X 10^^ C t7T~ 

t'oi l-0±^-2 l-9to-5 1.4tf^xl0-" l.stHxlO^" C T, 2 

45 2.0lS ;^ l-63tlJ l^^ 2.lt°oi x lO-^^ l.Ot"'^ x lO^^ X 3, 4 

5 50.0ti° l] 1.8^!^-2 1.7^0-8x10-12 5.1^^-^x1035 C T, 5 

3t°;J« 50.0;i°;;^ 0.78;°;°^ 2.9tlf x lO-i^ 8.6l3"J-4 x 103^ C T, 5 

EJ-^ 2.6t^3 4.0t2-e X 10-15 1.0^2-3 x 1030 C 6,7 

t|]-,i,2 1.8+0.2 2.lt}-[; X 10-13 1.6;t-ixl03" C T, 8 

to-03 2.lt^l 9.6tf^x 10-14 9.5™ X 1029 X 9,10 

2.5iO-o 2.6i^-J X 10-15 3.6^1^x10'' X 11,7 

+0-0^ 1.3to-2 4.0tfo X 10-12 4.2t5:|xl03i C T, 12 



1, 


q+0.1 
■^-0.1 


5. 


1+1.3 

-■--1.0 


X 


10- 


■14 


4, 


,2+1-3 X 1020 


X 


11, 13 


1, 


7+0.1 
■ ' -0.1 


7. 


0+2.8 
■^-1.8 


X 


10- 


■14 


4, 


.5tli X 1030 


X 


9, 14 


1, 


co+O.lO 
■"^-0.10 


6. 


n+2.0 

■«-1.5 


X 


10- 


■13 


3. 


,1+2-2 ^ ig33 


c 


T, 15 


1, 


•"-0.3 


6. 


c;+7.3 
^-3.3 


X 


10- 


-14 


5, 


,3+iV X 1031 


c 


16, 7, 14 


1, 


■^^-0.09 


3. 


4+1.3 
^-1.0 


X 


10- 


-12 


4, 


,5+0-^ X 1033 


c 


T, 17 



B0950+08 0.253 2.298 x lO'i^^ 0.03tg-^^ 0.262iO-005 
B1055-52 0.197 5.833 x lO'i^ 0.027+!^-^i^2 OJIqI 
J1124-5916 0.135 7.471 x 10-^3 0.28t[;-J^4 6.2t^-^ 
J1509-5850 0.089 9.170 x lO-i^ 0.8+1^-2 2.6li-i^ 
J1617-5055 0.069 1.351 x 10-i3 3.3;!^-^ '-^-^-is 

10-^ 2.0^0-^ 1.1^3-0x10-13 9.7t^4/xl03i C 18,19 

tli 1.4^01 3.313-3x10-12 9.9^342' X 1033 C T, 20 

til l.6t°oi 6.9^0-^x10-13 2.2t«-4xl033 C 21,7 

tli l.2t°oi 2.5tiVx 10-14 3.6+3^43 x 1031 C 22,7 

'^2:5 l-Otas 2.5t?;i X 10-12 7.5123^4 x 1033 c 23,24 

B1823-13 0.101 7.506 x 10-i4 l.Otal ^-^-le l-7to7 5.4tf.3 x 10-i4 9.8tg320-3 x 103i C T, 7 

J1846-0258 0.326 7.083 x IO-12 3.694' G-Otfo l-9ltoio '^■^-38 x 10^^^ S.itl^ x 1033 C T, 25 

B1853+01 0.267 2.084 x 10-^3 59 2.6;'}'i^ l-3to5 7.7^5^x10-14 6.2^337^x1031 C 26,27 

B1929+10 0.227 1.157 x lO-i^ 0.24l:0-io 0.361+!^-HJ2 3.01^'^ 5.3tS-5 x 10-i4 8.2^0:5 x 1020 C T, 2 



B1706- 


-44 


0. 


,102 


9, 


,298 


X 


10- 


■14 


0, 


.55 


2, 


,7: 


J1747- 


-2958 


0. 


,099 


6, 


.136 


X 


10- 


■14 


2, 


.615 


5, 


.0: 


B1757- 


-24 


0. 


,125 


1, 


.279 


x 


10- 


13 


3 


rn+1.30 
■^"-1.10 


5, 


.2: 


J1809- 


1917 


0. 


,083 


2, 


.554 


X 


10- 


■14 


0, 


.72 


3, 


.5: 


J1811- 


-1925 


0. 


,065 


4, 


.400 


X 


10- 


■14 


2, 


9+0.8 
■^-0.6 


5, 


.0: 



'-0.09 ""J^-^-COOS •-'•"-0.3 ^■'^-2.7 ^ °-^-4.3 

J1930+1852 0.137 7.506 x 10-i3 1.6 5.0^2^ l-2to;2 l-8t?i x IO-12 5.412^' x 1033 C T, 28 



B1951+32 0.040 5.845 x 10-i5 0.299 2.01^-^ l-64too9 2.0t{]1 x IO-12 9.7tff x 1032 C T, 29 
J2229+6114 0.052 7.827 x 10-i4 0.27tg;0^ O-SIm J-OSto'lo 3.7^1;^ x 10-i3 2.9t^-0 x 103i C T, 30 



Note. — C: Chandra, X: XMM-Newton 
The error bars of /x arc derived from the errors in the spectral index and the normahzation of spectrum fitting, the error bars of .Lxjpsr are 
derived from the errors of both the fluxes and the distances. If more than one papers is listed in Ref. column, the first one is -where the spectral 

properties arc from and the other ones are those where the distance is from. 

References. — T; This paper; [1] Roberts et al.(1993); [2] Chatterjee et al.(2004); [3] Willingale et al.(2001); [4] Trimble & Woltjer(1971); [5] 
Bradlcy(2007); [6] Ogelman & Tepedelenlioglu(2004); [7] Cordes & Lazio(2002); [8] Faherty et al.(2007); [9] De Luca et al.{2005); [10] Brisken et 
al.{2003); [11] Becker et al.(2004); [12] Caraveo et al.(2001); [13] Brisken et al.(2002); [14] Kramer et al.(2003); [15] Gaensler & Wallacc(2003); 
[16] Hui & Becker(2007); [17] Paron et al.(2006); [18] Gotthelf et al.(2002); [19] Koribalski et al.(1995); [20] Gaensler et al.(2004); [21] Kaspi et 
al.(2001); [22] Kargaltsev & Pavlov(2007) ; [23] Roberts et al.(2003); [24] Green et al.(1988); [25] Leahy & Tian(2007); [26] Petre et al.(2002); 
[27] Cox et al. (1999); [28] Lu et al.(2002); [29] Strom & Stappers(2000); [30] Kothes et al.(2006). 
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Table 2 

Emission properties of 24 PWNe 



PSR Name 


1022cm-2 


kpc 


r 


J0205+6449 




o 9+1.6 


2. 


B0355+54 


0.21°;? 


1 o+°-2 

^-•"-0.2 


1. 


B0531+21 


0.345 


2.0l°j 


2. 


J0537-6910 
B0540-69 


"•^"J~0.03 

0.43 


so.otin 


2. 
1. 



pwn 



fx (2-lOkeV) 

,pwn 



Det. Ref. 



"'-0.02 
9+0.3 

0.006 

43+0.06 
96^ 



ergs s cm 



ergss 



J06334 


-1746 





001 


B0833- 


-45 





017 


J1016- 


-5857* 





q+0.7 

"^-0.3 


J1124- 


5916 





281 


J1509- 


-5850 





O+0.9 
^-0.4 



J1617-5055 

B1706-44 

J1747-2958 

B1757-24 

J1809-1917 



O q + 0.3 

U.D_Q 2 
9 f^+Q.l 

^•o-o.i 

3.5 

0.72 



J1930+1852 
B1951+32 
J2021+3651* 
J2229+6114 



1.64 



+0.10 
-0.09 

o.aotHl 



0.27 



-0.1 
+0.08 
-0.07 



"•^■J-o.oe 
29+° °^ 

"•^^-0.05 

q 0+0.6 

-0.6 
0-^-0.9 

2.6+1° 



-1.0 



+ 1.6 
1.6 
+0.9 

-2.5 

5.2t|l 
3.5tl;l 



3.3 
2.7 
5.0 



J1811-1925 


2 


j^+O.l 


5 


B1823-13 


1 


n+8:6 

"-0.6 


3 


J1846-0258 


3 


7+0.1 
'-0.1 


6 


B1853+01 


5 





2 


B1929-M0 





24+0.10 
^^-0.09 






0.361 



+2.5 

+i:i 

-1.6 
+3.0 
-3.0 
+ 1.3 
-1.3 

+0.010 
-0.008 



+2.5 
2.5 



5.0 

8.0+^-° 



0.8 



-4.0 
+0.4 
-0.4 



-0.06 
h0.08 
-0.08 



q+0.5 

•-"-o.s 

C9+0.02 
"3^-0.02 
+0.7 
0.5 
+0.10 
0.10 
+0.8 
0.4 



9^ 

86"' 

3^ 



+0.3 
-0.2 
[+0.2 
'-0.3 
+0.10 
0.10 
0+0.6 

4^ 



93 



-0.6 
1+0.1 
-0.1 



7+0.2 
-0.2 

q + 0.4 

"^-0.5 

oq+0.05 

9+0.2 
^-0.2 
4+0.6 
*-0.5 



89: 



+0.10 
0.09 
74+0.03 
'*-0.03 
7+0.3 



-0.2 
+0.1 
-0.1 



5 


f.+0.3 
O-0.3 


X 


10- 


12 


6 


0-I-9.4 
^-5.2 


X 1033 


C 


T 


1 


4+1-2 
^-0.7 


X 


10- 


13 


1 


0+3.1 
^-1.2 


X 1031 


C 


T 


2 


1+0.1 
■■■-0.1 


X 


10- 


08 


1 


0+0.7 
"-0.5 


X 103^ 


X 


1 


2 


O+0.4 
^-0.4 


X 


10- 


12 


8 


4+5.6 
^-3.8 


X 1035 


c 


T 


5 


4+1.4 
^-1.1 


X 


10- 


12 


1 


fj + 1.3 
"-0.8 


X 1036 


c 


T 



9 q+13.4 

-1-8 

6.8ti7 
1 c+8.0 
1-<J-1.5 



X 10- 


14 


3 


8 


X 10- 


11 


6 


4 


X 10 


-13 


2 


5 


X 10- 


12 


3 


1 


X 10- 


13 


1 


2 



1 4+5-3 
1-^-1.0 
9 9+0.4 
^•^-0.3 

2 6+1'^ 

^•"-0.9 
Q 0+9.6 
"J-O-2.9 



X 10 

X 10- 
X 10" 
X 10 
X 10- 



-12 
-13 
-11 
-13 
-14 



+ 19.2 



X 10 



29 



tf2« X 1032 

til X 1034 
+ 1%" X 1032 



3 


7+2.7 
' -2.3 


X 


10- 


13 


4 


8 


4 


0+5.7 
"-2.2 


X 


10- 


13 


3 


5 


4 


f^+1.6 
0-1.2 


X 


10- 


12 


1 


4 


1 


q+3.5 
"^-0.9 


X 


10- 


13 


4 


1 


2 


q+0.8 
"^-0.6 


X 


10- 


13 


3 


4 



+ 13.8 

-4.4 
+ 11.6 
-2.8 
+ 2.8 



X 1032 
X 1032 
1.4:^^-^ X 1034 
tfs^ X 1032 
+5-j X 1032 

l.Ot^-^ X 1034 
2.5^2137 X 1032 
9.3+15/ X 1034 



-7.4 
1+5.1 



-132 



5 


4+1.6 
^-1.3 


X 


10- 


12 


1.6 


5 


1+0.4 

■■■-0.4 


X 


10- 


12 


2.5 


1 


1+0.5 
■■■-0.4 


X 


10- 


12 


8.4 


4 


0+2.0 
"-1.3 


X 


10- 


13 


3.0 



2.11°-^ X 103 
5.9tiy X 10^9 

tfl X 1034 

+3.5 X 1033 



-1 

+ 19.2 
-7.2 
+7.3 
-2.5 



X 10 
X 10 



33 



c 
c 
c 



31 



c 
c 
c 
c 
c 

c 
c 
c 
c 
c 

c 
c 
c 
c 



T 
T 

T, 2 



C T 
C 3 



T 
4 
T 
5 
6 

7 
T 
T 



T 
T 
9 
T 



Note. — C: Chandra, X: XMM-Ncwton 
★ PSRs .11016-5857 and J2021+3651 arc not listed in Tabic 1. Their periods are P = 0.107 and 0.104 s, respectively, and their period 

derivatives are P = 8.08 X 10-" and 9.56 X 10-"ss"\ respectively. 

The error bars of /x are derived from the errors in the spectral index and the normalization of spectrum fitting. The error bars of ix,pwn aJ^e 

derived from the errors of both the fluxes and the distances. 

References. — T: This paper; [1] Willingale et al.(2001); [2]Ruiz & May(1986); [3] Hui & Becker(2007); [4] Gotthelf et al.(2002); [5] Kaspi 
et al.(2001); [6] Kargaltsev & Pavlov(2007) ; [7] Roberts et al.(2003); [8] Petre et al.(2002); [9] Hessels et al.(2004). 
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Table 3 

Correlation coefficients between parameters 



y 


X 


r 




Ps 








bnc 


-£/X,psr 


P 


-0.73 


-0.66 


0.0002 


27.1±0.2 


-5.2±0.2 


28.9±0.7 


-3.6±0.7 




P 


0.70 


0.69 


0.0001 


49.7±0.8 


1.29±0.06 


51.0±3.9 


1.4±0.3 




B 


0.41 


0.46 


0.0165 


1.6±1.6 


2.5±0.1 


10.9±9.4 


1.7±0.8 




T 


-0.84 


-0.81 


< 10-4 


38.1±0.3 


-1.19±0.05 


38.9±0.9 


-1.4±0.2 




E 


0.89 


0.82 


< 10-4 


-0.8±1.3 


0.92±0.04 


-4.2±3.7 


l.OiO.l 


r 

psr 


P 


0.57 


0.60 


0.0010 


2.03±0.08 


0.49±0.08 


2.3±0.2 


0.7±0.2 




P 


-0.45 


-0.42 


0.0274 


0.3±0.3 


-0.09±0.02 


-1.5±1.3 


-0.23±0.09 




B 


-0.22 


-0.22 


0.2598 


2.5±0.6 


-0.08±0.05 


4.6±2.6 


-0.2±0.2 




T 


0.58 


0.44 


0.0205 


l.OiO.l 


0.12±0.02 


0.5±0.3 


0.26±0.07 




E 


-0.64 


-0.51 


0.0067 


5.7±0.6 


-0.11±0.02 


8.5±1.6 


-0.19±0.04 




C 


-0.64 


-0.50 


0.0081 


2.9±0.2 


-0.54±0.08 


3.8±0.5 


-0.9±0.2 




■t'X,psr 


-0.55 


-0.54 


0.0033 


5.4±0.5 


-0.12±0.02 


6.2±1.4 


-0.14±0.04 


-^X,pwn 


P 


-0.54 


-0.40 


0.0519 


29.8±0.3 


-3.8±0.3 


29.9±1.2 


-3.2±1.1 




P 


0.66 


0.61 


0.0015 


54.3±1.7 


1.6±0.1 


52.3±4.7 


1.5±0.4 




B 


0.43 


0.52 


0.0090 


21.2±3.0 


1.0±0.2 


11.5±9.8 


1.7±0.8 




T 


-0.82 


-0.76 


< 10-4 


42.4±0.5 


-2.1±0.1 


40.5±1.1 


-1.7±0.3 




E 


0.86 


0.75 


< 10-4 


-19.6±3.0 


1.45±0.08 


-14.9±6.0 


1.3±0.2 






0.94 


0.91 


< 10-4 


-1.9±3.2 


1.08±0.10 


-0.3±2.6 


1.02±0.08 


F 

^ pwn 


P 


-0.25 


-0.21 


0.3216 


1.14±0.05 


-0.65±0.03 


1.2±0.3 


-0.4±0.3 




P 


0.38 


0.40 


0.0551 


4.8±0.2 


0.22±0.01 


4.1±1.3 


0.2±0.1 




B 


0.27 


0.35 


0.0889 


-0.2±0.3 


0.18±0.03 


-1.5±2.4 


0.2±0.2 




T 


-0.45 


-0.45 


0.0279 


2.94±0.04 


-0.27±0.01 


2.5±0.4 


-0.22±0.09 




E 


0.45 


0.43 


0.0367 


-5.7±0.3 


0.203±0.008 


-4.3±2.4 


0.16±0.07 




r 

^ psr 


-0.01 


0.04 


0.8554 


0.2±0.3 


1.1±0.2 


1.6±0.3 


-0.01±0.20 




■t'X.pwn 


0.70 


0.70 


0.0001 


-4.9±0.7 


0.20±0.02 


-3.8±1.2 


0.16±0.04 






0.73 


0.75 


< 10-4 


3.2±0.2 


0.46±0.08 


2.6±0.2 


0.31±0.06 


/x,pwn/ /x,psr 


E 


0.13 


0.15 


0.4917 


-2.3±5.7 


0.07±0.15 


-20.6±43.4 


0.7±1.2 



Note. — r is Pearson correlation coefBcient, while r's and ps arc Spearman rank eorrelation coefficient and significance level, respectively. 
Coefficients ae and be (ftne and fone) are obtained in the linear fitting with (without) observational errors. 
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28 30 32 34 36 38 28 30 32 34 36 38 

Ig U.psr (ergs s'^) Ig U.pwn (ergs s"'') 

Fig. 1. — The relations between the non-thermal X-ray luminosity in 2-10 keV and the photon index for RPPs (left panel) and PWNe 
(right panel). The solid lines are the best LSM fit without observational errors taJjen into account, while the dotted line with observational 
errors. The dashed lines mark the cases of r=2 for comparisons. 




32 34 36 38 34 35 36 37 38 

Ig Edot (ergs s"'') Ig Edot (ergs s"'') 

Fig. 2. — Left: The non-thermal X-ray luminosities (ix,psr) in 2-10 keV from the Chandra and XMM-JVewton observations vs. spin-down 
power E of 27 pulsars. Right: The non-thermal X-ray luminosity (i^x.pwn) in 2-10 keV from the Chandra and XMM-Newton observations vs. 
spin-down power E of 24 PWNe. The bold solid line is the best fit to the data by LSM without observational errors included, while the thin 
solid line is the fit with observational errors. For comparison the resulted relations in previous works are also marked: the dash-dot-dot-dot 
line corresponds to results of Becker & Triimper (1997); the dash-dot line to Saito (1998); the short dash line to Possenti et al. (2002); and 
long dash line to Cheng et al. (2004). 
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-1.5 -1.0 -0.5 0.0 -15 -14 -13 -12 -11 3 4 5 6 7 
Ig P(s) Ig Pdot(s s"') Ig t (yr) 

Fig. 3. — Relations of Lx.psr vs. P, P and r. The solid lines are the best LSM fit without observational errors considered while the dotted 
lines are with errors. The fitting results axe also listed in Table 3. 




Fig. 4. — Relations of Fpgr vs. P, P, r and E. The solid lines are the best LSM fit without observational errors considered while the dotted 
lines are with errors. The fitting results axe also listed in Table 3. 
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Fig. 6. — Relations of Fpwn vs. P, P, r and E. The solid lines are the best LSM fit with observational errors considered while the dotted 
lines are without errors. The fitting results are also listed in Table 3. 
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■ I I I I I I ■ . 

-7 -6 -5 -4 -5 -2 -1 

U.pwn/Edot 

Fig. 7. — The relation between Fpwii and ix,pwn/^- The solid line is the best LSM fit without observational errors considered while the 
dotted line is with errors. The fitting results axe also listed in Table 3. 




30 32 34 36 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Ig U.psr (ergs s'^) Tpsr 

Fk;. 8. Left: the relation of Lx.psr vs. ix,pwn. The solid lines are the best LSM fit with the observational errors considered while 
the dotted lines are without errors. The fitting results are also listed in Table 3. The dash line marking the case of ix.psr = l^x,pvjn is for 
comparison. Right: the relation of Fpsr vs. Fpwn. 
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33 34 35 36 37 38 39 
Ig Edot (ergs s"'') 

Fig. 9. — The relation of /x,pwn//x,psr vs E. The dash line indicates /x,pwn = /x, 



